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ABSTRACT 



We have detected C91o; (8.5891 GHz) emission toward 4 ultra-compact Hll 
regions (UCHs; W49G, J, L & C) in the W49 North massive star forming region 
with the Very Large Array (VLA) at 3" resolution. No carbon line emission 
was detected toward UCHs W49F, A, 0, S and Q at this frequency to a 3a 
level of 2 mJy We also observed the same region in the C75a line (15.3 GHz) 
with no detection at a 3a level of 6 mJy with a 1".7 beam. Detection of line 
emission toward these sources add supporting data to the earlier result of Roshi 
et al. (2005a) that many UCHs have an associated photo-dissociation region 
(PDR). Similarity of the LSR velocities of carbon recombination lines and H 2 CO 
absorption toward UCHs in W49 North suggests that the PDRs reside in the 
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dense interface zone surrounding these Hn regions. Combining the observed 
carbon line parameters at 8.6 GHz with the upper limits on line emission at 15.3 
GHz, we obtain constraints on the physical properties of the PDRs associated 
with W49G and J. The upper limit on the number density of hydrogen molecule 
obtained from carbon line models is ~ 5 x 10 6 cm~ 3 . 

Subject headings: ISM: HII regions — ISM: general - radio continuum: ISM - 
radio lines: ISM - line: formation 



1. Introduction 

Observations of molecular material in the vicinity of ultra-compact H II regions (UCHs) 
have shown that the ambient molecular gas may have very high densities. For example, 
observations of high density molecular tracers in the W49 North molecular cloud core show 
that the molecular density is a few times 10 6 cm -3 (Serabyn, Glisten & Schulz 1993; see also 
Churchwell 2002 and references therein). The presence of such high density gas surrounding 
the UCHs has other effects. In particular, the far-ultra violet (FUV; 6 - 13.6 eV) photons 
that escape the UCH produce photo-dissociation regions (PDRs) at the interface between 
the ionized and molecular material ( Hollenbach & Tielens 1997). In these PDRs, carbon 
will be ionized since it has a lower ionization potential (11.3 eV). The physical conditions in 
such PDRs are ideal to produce observable carbon radio recombination lines (RRLs). 

A recent search for carbon RRLs toward UCHs have detected lines in a majority of 
directions (Roshi et al. 2005a). These observations were done near 9 GHz toward 17 UCHs 
and carbon RRLs were detected from 11 sources (65 % detection rate). The large detection 
rate indicates that most UCHs have an associated PDR. High-angular resolution, multi- 
frequency carbon RL observation from these PDRs are useful in addressing important issues 
related to UCHs and their evolution. For example, such a data set can be used : (1) to 
estimate the physical properties of the PDR (Roshi et al. 2005a, b; Garay et al. 1998, 
Natta, Walmsley & Teilens 1994) and (2) to address the 'lifetime problem' of UCHs (Wood 
& Churchwell 1989) by investigating whether some UCHs are pressure confined (De Pree, 
Rodriguez & Goss 1995; Roshi et al. 2005a). 

In this paper, we report the detection of C91« (8.5891 GHz) emission toward several 
UCHs in W49 North. The W49 star forming region lies close to the Galactic plane (/ = 
43°. 2, b = +0°.00) and is embedded in a giant molecular cloud of size ~ 100 pc (Simon et 
al. 2001), at a distance of 11.4±1.2 kpc (Gwinn, Moran & Reid 1992). Radio continuum 
observations have identified 45 UCHs within the molecular cloud (De Pree, Mehringher Goss 
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1997). The properties of these UCHs were well studied through radio continuum as well as 
hydrogen and helium RL observations (De Pree et al. 1997; De Pree et al. 2000). Here we 
present new carbon RL observations of W49 North with the Very Large Array (VLA). The 
details of the observations are given in §2. In §3, we provide the parameters of the observed 
carbon lines and present models for the line formation. A discussion of the results obtained 
is given in §4. 



2. Observation and data reduction 

The observations were carried out on 22 November 1998 using the Very Large Array 
of the National Radio Astronomy Observatory 1 . We used the C array for the observations. 
The helium and carbon RRLs at both 8.6 GHz ( He and C 91a ) and 15.3 GHz (He and C 
75a) were interleaved during the 8 hour period. The synthesized beam is 1".8 x 1".6 (PA: 
-18°) at 15.3 GHz and 3".l x 2". 9 (PA: -17°) at 8.6 GHz. The observations were centered 
at 25 km s _1 (LSR) to center the recombination lines between the He and C lines ( the 
velocity difference between He and C is about 30 km s _1 ). The observed field was centered 
at the coordinates RA (J2000) = 19M0 m 13 s .41 and S (J2000) = 09°0l'l6".l. We used a 
total bandwidth of 3.125 MHz with 127 channels spaced by 24.4 kHz (velocity resolution ~ 
1 km s _1 ) at 8.6 GHz and 6.25 MHz with 63 channels spaced by 97.7 kHz (velocity resolution 
~ 2 km s" 1 ) at 15.3 GHz. Two orthogonal circular polarizations were observed. The source 
3C286 was used as the absolute flux density calibration and phase calibration was carried 
out by frequent observations of the compact source B1920+154. The bandpass response was 
determined by observations of 3C84 and 3C273. 

The images were made using natural weighting in order to improve the sensitivity to 
extended sources. The continuum was subtracted using UVLSF in AIPS based on the line 
free channels. The final line cubes have an RMS noise of 0.7 mJy /beam (Hanning smoothed) 
at 8.6 GHz and 1.9 mJy /beam at 15.3 GHz. 



1 The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under a cooperative agreement by Associated Universities, Inc. 
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3. Results 

3.1. Carbon recombination line emission 

The central image in Fig 1 shows the 8.6 GHz continuum at ~3" resolution. The RMS 
noise in the continuum image is 7 mJy/beam. The spectra at 8.6 GHz obtained toward 
six sources are also shown in figure 1. The C91a line was detected toward four UCHs 
(W49G, J, C & L). The central coordinates and the extent of the region over which the 
data are averaged are listed in Table 1. Table 2 gives the parameters of both carbon and 
helium lines obtained from the six spectra. The strongest carbon line (flux density 6 mJy) is 
observed toward W49G. This result is expected since the carbon line emission is dominated 
by stimulated emission at these frequencies (Natta et al. 1994; Roshi et al. 2005a, b) and 
W49G is the strongest continuum source in the observed field. Weak carbon line emission is 
observed toward W49J, C and L. The parameters of carbon line emission toward W49L may 
be affected by the larger velocity extent of helium line emission. No carbon line emission is 
detected in our 15.3 GHz data to a 3cr level of 6 mJy with a beam size of I". 7. 



3.2. Comparison with molecular line emission 

We can compare carbon RRL emission with existing CS (Serabyn et al. 1993) and H 2 CO 
(Dickel & Goss 1990) observations. Multi-transitional observations of CS have provided 
information on the structure of molecular cloud in W49 North and its kinematics. The cloud 
core consists primarily of three molecular clumps, designated as CS-NE, CS-C and CS- 
SW, with sizes between 0.4 - 1 pc. Combining the CS observations with H2CO absorption 
studies, Serabyn et al. (1993) suggested that W49G, C and F are associated with the 
clump CS-C; W49A and L are probably associated with CS-SW and CS-NE, respectively. 
The association of W49J with any one of three molecular clumps is not certain. Here 
we focus on the associations of UCHs in W49 North which are relevant for the present 
discussion. The LSR velocity of CS lines from clump CS-C is 4.3±1.0 km s -1 , similar 
to the velocity of the carbon RRLs (4.2±1.0 km s _1 ) in the direction of W49G and C. 
These two sources also exhibit H 2 CO absorption components with an average velocity of 
5.2±0.8 km s _1 . The carbon RRLs observed from W49L and J have velocities of 7.1±0.6 
km s _1 , different from the velocity of the CS line from the clump CS-NE (11.9±1.0 km s _1 ). 
The H 2 CO absorption components observed toward these sources have velocities 7.2±0.2 
km s _1 (Dickel & Goss 1990; We estimated approximate errors on the central velocity of 
H 2 CO absorption line from the data given in the reference). The similarity of the velocity of 
the H 2 CO absorption components and carbon RRL emission toward all four sources suggests 
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that the PDR producing carbon lines reside in the gas responsible for H 2 CO absorption. 

3.3. Models for carbon line emission 

We present models for carbon line emission toward W49G and J. Constraints on the 
physical properties of the PDR could not be obtained toward W49C because the line emission 
is weak. The W49L data are excluded from modeling since the line parameters may be 
affected by the large velocity extent of the helium line emission. The detailed steps involved 
in modeling the carbon line emission are described by Roshi et al. 2005b. A homogeneous 
'slab' of PDR material placed in front of the UCH is considered for modeling. The radiative 
transfer equation for non-LTE cases is solved to obtain the RL flux density. We used the 
computer code originally developed by Brocklehurst & Salem (1977) and later modified by 
Walmsley & Watson (1982) and Payne, Anantharamaiah & Erickson (1994) to calculate the 
non-LTE departure coefficients. 

The abundance of carbon required for the calculation of departure coefficient is taken as 
3 x 10~ 4 , which implies a depletion factor of 25 % (Natta et al. 1994). The departure coef- 
ficients also depend on the background radiation field, assumed to be a thermal background 
from the UCH. The continuum emission at 8.6 and 15.3 GHz from the UCH (see Table 2) is 
used to estimate an average emission measure for an assumed temperature of 9000 K. The 
estimated emission measures for W49G and J are summarized in Table 3. For a given PDR 
gas temperature, we varied the PDR emission measure to obtain the carbon line flux densi- 
ties which are consistent with the 8.6 GHz detection and 15.3 GHz upper limit. Fig 2 shows 
the results of modeling for W49G and J. The data at the two frequencies provide an upper 
limit on the PDR density and lower limit on the line of sight thickness of the PDR. These 
limits are obtained for assumed temperatures of 500 and 1000 K (Natta et al. 1994, Roshi 
et al. 2005b) and summarized in Table 3. The electron density obtained from modeling is 
converted to molecular density n# 2 using the assumed abundance of gas phase carbon and 
by assuming that hydrogen is in molecular form in the PDR. 

4. Discussion 

The detection of carbon RRLs toward a number of sources in the star forming region 
W49 North confirm the conclusion (Roshi et al. 2005a) that most UCHs may well have an 
associated dense PDR. These PDRs are formed at the interface between the ionized and dense 
neutral material, perhaps a dusty zone separating the two region as in the case of OH maser 
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sources W49B and G (Dickel & Goss 1990). These interface zones are dense with H 2 densities 
~ 10 5 cm" 3 . Observational evidence for the presence of such zones surrounding W49C and 
J exists; for W49L the existence of a distinct interface region is uncertain. Dickel & Goss 
(1990) suggest that the interface zone is responsible for the H 2 CO absorption. Carbon 
RRLs observed toward the UCHs in W49 North have LSR velocities similar to that of 
H2CO absorption components. This similarity suggests that the PDR responsible for carbon 
line emission resides in these interface zones surrounding UCHs. Modeling the carbon line 
emission has provided upper limits on H 2 densities, which is ~ 5 x 10 6 cm -3 . This upper 
limit is consistent with the densities inferred for the molecular gas in the interface zone from 
H 2 CO absorption studies (Dickel & Goss 1990). 
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Table 1. Details of the region over which the spectra are averaged 



Source RA(J2000) a DEC(J2000) a Avg. area a 





1 


'h m s'j 


( 


o ' 


") 


(arcsec 2 ) 


W49G 


19 


10 


13.54 


09 


06 


12.6 


19 


W49J 


19 


10 


14.15 


09 


06 


15.2 


14 


W49L 


19 


10 


14.44 


09 


06 


21.2 


16 


W49C 


19 


10 


13.13 


09 


06 


18.6 


11 


W49F 


19 


10 


13.34 


09 


06 


22.1 


9 


W49A 


19 


10 


12.84 


09 


06 


12.1 


12 


W490 


19 


10 


16.35 


09 


06 


07.3 


18 


W49S 


19 


10 


11.74 


09 


05 


27.0 


27 


W49Q 


19 


10 


10.56 


09 


05 


13.9 


50 



a The spectrum for each source in Fig. 1 and line 
parameters in Table. 2 are obtained by averaging 
the data over the area given here. The central 
coordinates of the area are also given. 
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Table 2. Parameters of the observed recombination lines 



Source 


fobs 


S c a 


Line 




AV h 


Vlsr 




(GHz) 


(Jy) 




(mJy) 


(km s ) 


(km s ) 




o.D 


i.y^u.o j 


L/yio 


o.y 


y.o ^u.oj 


a iVn 
4.11^11.0 j 








tteyia 


D. < \y)-<j) 


lo.o^u.y j 


£. Aid A\ 
D.4(U.4 J 




10. o 


Z.4(U.O J 


O / OOi 


(P. A\ 






W4yj 


O.D 


U.o^U.l J 


L/yia 


o.l (U.zJ 


10.1(1.0 J 


7 1 /n (i\ 








T4oQ1 ru 

neyio 


(.1 (U.ZJ 


ly.o^u.o ) 


D.U^U.O ) 




10. 


U.O^U. 1 ) 


o / ou 


{6.0) 








O.D 


U.O^U.l J 


oy 10 


Q QC 9"\ 


ic 7/1 rA 

id. ( (i.yj 


7 ifn k\ 








He91a 


10.3(0.2) 


20.7(0.7) 


4.1(0.3) 


W49C 


8.6 


0.2(0.1) 


C91a 


l.l d (0.2) 


8.0 


4.3(0.9) 








He91a 


2.4(0.2) 


19.8(1.4) 


7.2(0.4) 


W49F 


8.6 


0.13(0.03) 


C91a 


(0.4) 












He91a 


2.9(0.1) 


16.2(0.9) 


7.1(0.4) 


W49A 


8.6 


0.2(0.1) 




(0.7) 






W490 


8.6 


0.5(0.1) 


C91a 


(0.7) 












He91a 


4.2(0.2) 


29.6(1.4) 


5.9(0.6) 


W49S 


8.6 


0.5(0.1) 


C91a 


(0.7) 












He91a 


2.2 d (0.2) 


26.4(2.6) 


-3.0(1.1) 


W49Q 


8.6 


0.5(0.1) 


C91a 


(1.3) 












He91a 


7.5(0.2) 


17.4(1.1) 


8.7(0.5) 



a Contimmm flux density with la error obtained over the area given in 
Table 1. 

b The line parameters are obtained from the spectrum averaged over the 
area given in Table 1. The errors on the parameters are la values. 

c Carbon line parameters are affected by the large velocity extent of the 
helium line emission. 

d Tentative detection. 
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Table 3. Results of modeling 



Source 


EM bg a 
(pc cm -6 ) 


TpDR 

(K) 


„PDR 
(cm -3 ) 


I 

(x 10~ 4 pc) 


n H2 
(x 10 6 cm- 3 ) 


W49G 


6.5xl0 7 


500 


< 3000 


> 0.7 


5.0 






1000 


< 3500 


> 1.9 


6.0 


W49J 


1.2xl0 7 


500 


<2800 


> 3.6 


4.5 






1000 


<3000 


> 10.1 


5.0 



a Estimated emission measure of the UCH using the continuum data 
for an assumed electron temperature of 9000 K. 
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Fig. 1. — 8.6 GHz continuum image of W49 North and spectra toward 6 UCHs. The angular 
resolution of the image is 3".l x 2". 9 (beam position angle —17°). The RMS noise in the 
continuum image is 7 mJy/beam. The contour levels in the image are (1, 2, 3, 4, 5, 6, 7, 8, 
9, 10, 11, 12, 13) x 100 mJy/beam and the grey scale ranges from —0.041 to 1.332 Jy/beam. 
We adopt the naming convention followed by Dickel & Goss (1990) to designate the UCHs. 
The spectra are obtained by averaging the data over a region specified in Table 1. The fitted 
Gaussian components to the line emission are also shown in each spectrum. Carbon line is 
detected toward W49G and weak lines are present in the spectra toward W49J, C and L. 
Spectra toward W49F and A are examples of our non-detections of carbon lines. 
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W49G W49J 




Frequency (GHz) Frequency (GHz) 



Fig. 2. — Carbon line flux density as a function of frequency from PDR models (see §3.3) 
toward W49G (left) and J (right). The solid and dashed curves correspond to gas temper- 
atures of 500 and 1000 K respectively. The electron densities and path lengths for these 
models are given in Table 3. The upper limits on carbon RL emission at 15.3 GHz and C91o: 
flux densities with ± la error are also marked. 



